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ABSTRACT 



Astrophysical gamma-ray spectroscopy is a most 
valuable tool for studying nuclear astrophysics, as 
well as recent star formation in the Milky Way. Af- 
ter a short, historical, introduction to the field, I 
present a brief review of the most important cur- 
rent issues. Emphasis is given to radioactivities pro- 
duced by massive stars and associated supernova ex- 
plosions, and in particular, those related to obser- 
vations presently carried out by INTEGRAL: short- 
lived ^-^Ti from CasA and SN1987A and long-lived 
^^Al and ^"Fe from massive stars; various candidate 
sources of positrons for the 511 keV emission of the 
Galactic bulge are also critically discussed. 

Key words: Nucleosynthesis, supernovae, radioactiv- 
ities, Milky Way, gamma-ray lines. 



1. HISTORICAL OVERVIEW 



The starting point of 7-ray line astronomy with cos- 
mic radioactivities is usually considered to be the 
landmark paper of Clayton, Colgate and Fishman 
(1969). That work clarified the implications of the 
production of ^®Ni (a doubly magic, and yet unstable 
nucleus) during explosive Si-burning in supernovae 
(SN). In particular, it opened exciting perspectives 
for 7-ray line astronomy, by suggesting that any su- 
pernova within the local group of galaxies woud be 
detectable in the characteristic 7-ray lines resulting 
from the radioactive decay of ^^Ni and its daughter 
nucleus ^^Co. 

In the 70's D. Clayton identified most of the ra- 
dionuclides of astrophysical interest (i.e. giving a 
detectable 7-ray line signal); for that purpose, he 
evaluated their average SN yields by assuming that 
the corresponding daughter stable nuclei are pro- 
duced in their solar system abundances. Amazingly 
enough (or naturally enough, depending on one's 
point of view) his predictions of average SN radionu- 
clide yields (Table 2 in Clayton 1982) are in excellent 
agreement with modern yield calculations, based on 



full stellar models and detailed nuclear physics (see 
Fig. 1). Only the importance of ^^Al escaped Clay- 
ton's (1982) attention, perhaps because its daughter 
nucleus ^^Mg is produced in its stable form, making 
the evaluation of the parent's yield quite uncertain. 
That uncertainty did not prevent Arnett (1977) and 
Ramaty and Lingenfelter (1977) from arguing (on 
the basis of Arnett's (1969) explosive nucleosynthe- 
sis calculations) that, even if only 10~^ of solar ^^Mg 
is produced as ^^Al, the resulting Galactic fiux from 
tens of thousands of supernovae (during the '-^l Myr 
lifetime of ^^Al) would be of the order of 10^'* cm~^ 



In the case of ^^Al nature appeared quite generous, 
providing a 7-ray fiux even larger than the opti- 
mistic estimates of Ramaty and Lingenfelter (1977): 
the HEAO-3 satellite detected the corresponding 1.8 
MeV line from the Galactic center direction at a 
level of 4 10"^ cm-^ s-i(Mahoney et al. 1984). 
That detection, the first ever of a cosmic radioac- 
tivity, showed that nucleosynthesis is still active in 
the Milky Way; however, the implied large amount 
of galactic ■^^Al (~2 M0 per Myr, assuming steady 
state) was difficult to accomodate in conventional 
models of galactic chemical evolution if SN were the 
main ^^Al source (Clayton 1984), since ^^Al would 
be overproduced in that case; however, if the "closed 
box model" assumption is dropped and infall is as- 
sumed in the chemical evolution model, that diffi- 
culty is removed, as subsequently shown by Clayton 
and Leising (1987). 

Another welcome mini-surprise came a few years 
later, when the ^^Co 7-ray lines were detected in 
the supernova SN1987A, a ~20 Mg star that ex- 
ploded in the Large Magellanic Cloud. On theoreti- 
cal grounds, it was expected that a SNIa (exploding 
white dwarf of ^1.4 Mq that produces ~0.7 Mq of 
^^Ni) would be the first to be detected in 7-ray lines; 
indeed, the large envelope mass of massive explod- 
ing stars (~10 Mq) allows only small amounts of 
7-rays to leak out of SNII, making the detectability 
of such objects problematic. Despite the intrinsically 
weak 7-ray line emissivity of SN1987A, the proxim- 
ity of LMC allowed the first detection of the tell- 
tale 7-ray line signature from the radioactive chain 
^^Ni^^^Co^^^Fe; this confirmed a 20-year old con- 
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Figure 1. Yields of radioactive nuclei from massive exploding stars. On the left, Clayton's (1982) predictions 
for average SN yields are given (based on the assumption that the corresponding stable daughter isotopes are 
produced in their solar abundances as radioactive progenitors); the uncertainty in ^^Fe yield stems from the 
unknown percentage of its contribution to the production of stable Ni (which may also be produced as unstable 
Co; the limits of the error bar correspond to percent and to 100 percent contribution, respectively). In the 
middle, recent theoretical results are plotted as a function of stellar mass (filled symbols, from Rauscher et al. 
2002), while open parallelograms indicate observational constraints from SN1987A, a 18-20 Mq star. On the 
right appear other observational constraints, either from CasA (for ^'^Ti) or from extragalactic SN (for ^^Ni, 
see also Fig. 3). 



jecture, namely that the abundant ^^Fe is produced 
in the form of radioactive ^^Ni. 

Those discoveries laid the observational foundations 
of the field of 7-ray line astronomy with radioactivi- 
ties. The next steps were made in the 90ies, thanks 
to the contributions of the Compton Gamma-Ray 
Observatory (CGRO). First, the OSSE instrument 
aboard CGRO detected the ^'^Co 7-ray lines from 
SN1987A (Kurfess et al. 1992); the determination of 
the abundance ratio of the isotopes with mass num- 
bers 56 and 57 offered a unique probe of the physi- 
cal conditions in the innermost layers of the super- 
nova, where those isotopes are synthesized (Clayton 
et al. 1992). On the other hand, the COMPTEL 
instrument mapped the Miky Way in the light of the 
1.8 MeV line and found irregular emission along the 
plane of the Milky Way and prominent "hot-spots" in 
directions approximately tangent to the spiral arms 
(Diehl et al. 1995); that map implies that massive 
stars (SNII and/or WR) are at the origin of galactic 
2^A1 (as suggested by Prantzos 1991, 1993) and not 
an old stellar population like novae or AGB stars. 

Furthermore, COMPTEL detected the 1.16 MeV 
line of radioactive ^^^j jj^ ^j^g Cas-A supernova rem- 
nant (lyudin et al 1994). That discovery offered 
another valuable estimate of the yield of a radioac- 
tive isotope produced in a massive star explosion (al- 
though, in that case the progenitor star mass is not 
known, contrary to the case of SN1987A). On the 



other hand, it also created some new problems, since 
current ID models of core collapse supernova do not 
seem able to account for the yield inferred from the 
observations; however, recent axisymmetric models 
of rotating star explosions offer interesting perspec- 
tives in that respect (see next section). 

Finally, after CGRO and before INTEGRAL, an- 
other important discovery was made in the field: the 
RHESSI experiment confirmed the detection of ^^Al 
in the Galaxy, but it also detected the character- 
istic decay lines of ^°Fe, another long-lived isotope 
with ~2 Myr period (Smith 2004, and this volume). 
The detected ^"Fe emissivity is, however, consider- 
ably lower than expected from current models of nu- 
cleosynthesis in core collapse SN. The implications 
of that discovery are studied in Prantzos (2004) and 
in Sec. 4. 



Table 1 summarizes our current picture of stellar 
radioactivities and associated 7-ray lines. In the 
following we shall focus on the radioactivities pro- 
duced by massive stars and associated supernova ex- 
plosions, and in particular, those related to observa- 
tions presently carried out by INTEGRAL. Similar 
reviews have been presented by Diehl and Timmes 
(1998) and Knodlseder and Vedrenne (2001). The 
current status of the physics of supernova explosions 
is reviewed by Hillebrandt (this volume). 
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Table 1. Important stellar radioactivities for gamma-ray line astronomy 



DECAY 
CHAIN 


MEAN LIFE* 


LINE ENERGIES (MeV) 
(Branching Ratios) 


SITE 
[Detected] 


NUCLEAR 
PROCESS 




0.21 


0.478 (0.1) 


Novae 


Expl.H 


56|\|j^56(-Q+^56pg 


0.31 


0.847 a.) 1.238 CO. 68) 
2.598 (0.17) 1.771 (0.15) 


SN 
[SN1987A] 
[SN1991T] 


NSE 




1.1 


0.122 CO. 861 0.136 CO. Ill 


SN 
[SN1987A] 


NSE 




3.8 


1.275 (1.) 


Novae 


Expl.H 




89 


1.157 ri.i 

0.068 (0.95) 0.078 (0.96) 


SN 

[CasA] 


a-NSE 


26A| + ^26Mg 


1.1 10*^ 


1.809 (1.1 


WR, AGB 
Novae 
SNII 
[Galaxy, Vela] 


St.H 
Expl.H 
St.Ne, Expl.Ne, u 


60pe^60Co^60|sjj 


2.2 10*5 


1.332 (1.) 1.173 (1.) 


SN 


n-capt 



[Galaxy] 



+ : positron emitters (associated 511 keV line) 

* : Double decay chains: the longest lifetime is given; Underlined : lines detected 
In parentheses : branching ratios; In brackets : sites of lines detected 
St.(ExpL) : Hydrostatic(Explosive) burning; NSE : Nuclear statistical equilibrium 
en: CK-rich "freeze-out"; n-capt : neutron captures; I^: neutrino-process 



2. 56^1 AND 44ti from CORE COLLAPSE 
SUPERNOVAE 

Both ^^Ni and ^''Ti are produced in the innermost 
layers of core cohapse SN, through explosive Si- 
burning. In the high temperatures resulting from 
the passage of the shock wave, matter enters Nuclear 
Statistical Equilibrium (NSE). Electron captures are 
too slow to change the neutron/proton ratio during 
the ~1 s timescale of explosive nucleosynthesis and 
thus most of the original composition (consisting of 
^'^Si, with equal numbers of protons Z and neutrons 
N) turns into ^^^Ni (Z=N=28). 

After the passage of the shock wave, material 
brought to NSE cools down. In environments with 
relatively low densities this happens in the presence 
of a large concentration of free a particles, which 
have no time to merge back to Fe-peak nuclei through 
the inefficient S-a reaction (the rate of which scales 
with density squared). This "a-rich freeze-out" pro- 
cess favours in particular the production of ■*^Ti(see 
discussion in Thielemann et al. 1996, and in partic- 
ular Figs. 3 and 4a). 

The yields of ^^Ni, ""^^Ti and other Fe-peak nuclei 
are extremely difhcult to evaluate from first princi- 



ples, at least in the framework of current models of 
core collapse supernova. The layers undergoing ex- 
plosive Si-burning are very close to the "mass-cut" , 
that fiducial surface separating the supernova ejecta 
from the material that falls back to the compact ob- 
ject (after the passage of the reverse shock). Since 
no consistent model of a core collapse supernova ex- 
plosion exists up to now (e.g. Janka et al. 2003), 
the position of the mass-cut is not well constrained. 
Current ID simulations suggest that the yield of *'*Ti 
is even more sensitive to the mass-cut than the one 
oi^^m (see Fig. 2). 

The presence of ^^Ni in SN1987A has been unam- 
biguously inferred from the detection of the 7-ray 
lines of the decay of its daughter nucleus ^^Co. The 
yield of ^®Ni has been estimated from the extrapo- 
lation of the early optical lightcurve to the origin of 
the explosion (precisely known thanks to the neu- 
trino signal, see Arnett et al. 1989 and references 
therein); the derived value, 0.07 M0, is often taken 
as a "canonical" one for core collapse SN. 

It turns out, however, that core collapse SN display 
a wide range of '^^Ni values, spanning a range of at 
least one order of magnitude, as can be seen in Fig. 
3 (Hamuy 2003 and references therein). Despite the 
large error bars, it seems that there is a clear cor- 
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Figure 2. ^^Ni and ^^Ti abundance profiles in the 
ejecta of a 20 Mq star (from models calculated by 
Rauscher et al. 2002). 



relation between the amount of ^^Ni and the energy 
of the explosion (obviously, because a shock of larger 
energy heats a larger amount of material to NSE con- 
ditions). 

'*'*Ti has not been directly detected in SN1987A up 
to now. Modelling of the late lightcurve of that 
supernova suggests that it may be powered by 1- 
2 10"'* Mq of '''^Ti (Motizuki and Kumagai 2003a) 
and similar values are obtained through analysis of 
the infrared emission line of the ejecta (Fransson and 
Kozma 2002). Note that the evaluation of the ""^Ti 
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Figure 3. Total ejected mass (upper panel) and ^^Ni 
yield (lower panel) as a function of ejecta kinetic en- 
ergy for core collapse supernovae (filled dots: SNII, 
crosses SNIb/c), according to Hamuy (2003). 



O 
E 
o 
o 

m 



1 


1 


1 1 1 

"Optimistic (Upper Limit)" Case 
<«Ti/=''Ni> = 2.9 
d = 50 kpc 


\^ 




- 


W. ESO 

\ ■» * 






\ '"f * * 
\ i\ Tj 

- \ 'Mr 1 




Collision to the Inner Ring 

? 

t„2.57y 


\r 

\\ ^ 


'-CTIO 


t t„2.63y 
• 


1 ?'Co 


1 


1 1 1 



Days after explosion 

Figure 4- Contributions of various radioactivities to 
the observed UVOIR lightcurve of SN1987A (from 
Motizuki and Kumagai 2003a). The late lightcurve 
requires the presence of ^2 10~'^ Mq of "^"^Ti; the 
derived '^'^Ti/'^Ni ratio is ^ 3 times higher than 
the solar ratio of the corresponding stable isotopes 
(^^Ca/>'^Fe)Q. 



yield through these methods suffers from consider- 
able uncertainties, due to the complex physics of 
the supernova heating and cooling and the role of 
positrons. The derived amount (and especially the 
higher value) is rather high compared to the results of 
ID models of 18-20 Mq stars (e.g. from Rauscher et 
al. 2002, see Fig. 1). Moreover, the '^'^Ti/^'^m mass 
ratio is around 3 times the solar ratio of the corre- 
sponding stable isotopes C^Ca/^^Fe)©, too high to 
be explained by current ID models, as will be dis- 
cussed below. 

Despite these high values, ''^Ti from SN1987A seems 
beyond the detection capabilities of INTEGRAL, 
now that the on-flight performance of SPI is well 
established. The analysis of Motizuki and Kumagai 
(2003a) indicates that the expected flux in the 1157 
keV line is ~5 10^® ph/cm^/s, i.e. considerably lower 
than the ^2 lO"'"^ ph/cm^/s sensitivity of SPIior an 
exposure of 1 Ms. 

The 7-ray lines of ''^Ti have been detected in the 
^^320 yr old Gas A supernova remnant, lying at a 
distance of ^3.4 kpc from earth. Both the high en- 
ergy line at 1.157 MeV and the low energy ones, at 
68 and 78 keV, have been detected, respectively by 
COMPTEL (lyudin et al. 1994) and Beppo-SAX 
(Vinketal. 2001). The detected flux of 3.3±0.6 10^^ 
ph/cm^/s from COMPTEL, points to a ^''Ti yield of 
~1.7 lO"'' Mq. Similar values, i.e. 1-2 lO"'' M©, 
are obtained through a study of the combined fluxes 
of the low energy lines (Vink and Laming 2002), al- 
though the modelisation of the underlying continuum 
spectrum makes the analysis very difficult. 

The '*'*Ti yield of CasA is comparable to the one 
implied by the optical and infrared observations of 
SN1987A. Note, however, that the CasA yield suf- 
fers from uncertainties related to the ionisation stage 
of the GasA remnant. ^''Ti decays by orbital elec- 
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Figure 5. Yield of^^Ti vs yield of^^Ni, from models 
and observations. Model results are from Limongi 
and Chieffi (2003, filled dots, with large variations 
in yields due to variations in both stellar mass - 
from 15 to 35 Mq - and explosion energy), Rauscher 
et al. (2002, crosses, for stars in the 15 to 25 
Mq range and explosion energies of llfi^ ergs) and 
Maeda and Nomoto ( 2003, asterisks, for axisymmet- 
ric explosions in 25 and 40 Mq stars, producing high 
^'^Ti/'^Ni ratios). '^'^Ti detected in CasA appears as 
a horizontal shaded band (assuming that its decay 
rate has not been affected by ionisation in the CasA 
remnant, otherwise its abundance should be lower, 
according to Motizuki et al. 1999). The amount 
of ^^Ti in SN1987A is derived from its late optical 
lightcurve (Motizuki and Kumagai 2003, see Fig. 4)- 
The diagonal dotted line indicates the solar ratio of 
the corresponding stable isotopes Ca/ Fe)Q. 



tron capture and an ionised medium could slow down 
its decay (Mochizuki et al. 1999); the derived yield 
could be considerably smaller in that case (Motizuki 
and Kumagai 2003b). On the other hand, ionisa- 
tion is expected to play a small role in the aformen- 
tioned derivation of the '^''Ti yield in SN1987A, but 
the remnant of that SN is rapidly evolving (Michael 
et al. 2002); if ^"'Ti enters a high ionisation stage 
in the near future, the expected 7-ray line fluxes 
from SN1987A should be accordingly reduced (see 
Motizuki and Kumagai 2003b). 




Figure 6. Abundance profiles of major isotopes as a 
function of ejecta velocity, after the bipolar explosion 
of a rotating star of 4O Mq; profiles have been aver- 
aged over all angles (from Maeda and Nomoto 2003). 



a more limited range of stellar masses (from 15 to 
25 Mq), but they considered only a fixed "canon- 
ical" value for the kinetic energy of the explosion; 
their results supersede those of Woosley and Weaver 
(1995), obtained with the same stellar evolution code 
but with an older set of nuclear reaction rates. Both 
calculations of RHHW02 and LC03 are made with 
the same reaction rate libraries and concern stars of 
solar initial metallicity. It is questionable to what ex- 
tent they apply to SN1987A, the progenitor of which 
presumably had an LMC metallicity of 0.3 Zq. 

The results of those calculations are plotted as ^''Ti 
yield vs ^^Ni yield in Fig. 5, where the solar ratio of 
the corresponding stable isotopes is also displayed as 
a diagonal line. It can be seen that: 

- There is broad agreement between LC03 and 
RHHW02, at least around the value of 0.07-0.1 Mq 
of the 56Ni yield. 

- Except for the lowest ^^Ni yields, '*''Ti and ^^Ni 
yields of LC03 both increase with the energy of the 
explosion, at a quasi-constant average ratio of ~3 
10-^ 



In summary: from optical observations we have a 
wide range of values for the ^^Ni yields of core 
collapse SN, and a precise value of 0.07 Mq for 
SN1987A; and for ''''Ti yields we have similar val- 
ues, i.e. 1-2 10"^ Mq, for both SN1987A (indirectly, 
through the modelisation of the UVOIR light) and 
for CasA (directly, through 7 ray lines, albeit with 
a systematic uncertainty resulting from poorly con- 
strained ionisation effects). How do these results 
compare with theoretical expectations? 

The most detailed relevant calculations have been re- 
cently performed for ID models, concerning stars in 
the 15-35 Mq range and parametrised by the energy 
of the explosion (Limongi and Chiefh 2003, LC03). 
Also, Rauscher et al. (2002, RHHW02) investigated 



- This value is ~3 times lower than the solar ratio 
of ('^'*Ca/56Fe)Q ^lO'^. This implies that such ex- 
plosions cannot produce the solar "^^Ca, since ^^Fe 
would be overproduced in that case (e.g. Timmes 
et al. 1996). Moreover, there is another important 
source of Fe, SNIa, which produce about 0.5-0.65 of 
solar ^^Fe, but very little '^''Ca; this makes the deffi- 
ciency of ''^Ca from core collapse SN even more se- 
rious than appearing in Fig. 5, since it implies that 
core collapse explosions should produce a "''^Ti/^^Ni 
ratio at least twice solar in order to compensate for 
the ^^Fe production of SNIa. Such a high ratio is 
also what is required to explain the high value of the 
44Ti yield in SN1987A. 

- In the case of CasA, one has two choices. The first is 
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to assume a very energetic explosion and, according 
to ID models, get a ^'*Ti yield compatible with CasA 
observations; but then a high ^^Ni yield of >0.2 M0 
is implied by those same models, leading to the ques- 
tion "why such a bright explosion as CasA went un- 
detected back in the 1680s?" (e.g. Hartmann et al. 
1996). The second option is to assume that the ^^Ni 
yield was much smaller than suggested by ID models 
producing high '^'^Ti (thus explaining why CasA went 
undetected); in that case it becomes clear that cur- 
rent ID models are unable to explain the situation 
in either CasA or SN187A. 

A solution to this problem of "missing ^■^Ti (and 
''^Ca)" is provided by multi-dimensional models of 
the energetic explosions of rotating massive stars 
(hypernovae). In such explosions, material falling 
onto the central remnant forms an accretion disk, 
and a fraction of that material is ejected in the 
form of two jets (or a bipolar wind), which inter- 
act with the stellar mantle. According to Nagataki 
et al. (1998), material along the jet axis undergoes 
higher temperatures and entropies (i.e. lower den- 
sities) than material in normal spherical explosions. 
Such conditions are particularly favourable to an "a- 
rich freeze-out" and lead to the production of large 
^^Ti amounts and '*'*Ti/^^Ni ratios. In Fig. 5 it can 
be seen that recent nucleosynthesis results obtained 
in the framework of such models (from Maeda and 
Nomoto 2003) may cure all the aforementioned prob- 
lems. 

At least in the case of SN1987A, there is convincing 
evidence that the explosion was axisymmetric and 
not spherical (see Wang et al. 2002 and references 
therein). On the other hand, ^''Ti produced in the 
jets is ejected at much higher velocities than in spher- 
ical symmetric explosions, reaching up to 25000 km/s 
according to Maeda and Nomoto (2003 and Fig. 6). 
This makes its lines rather broad and their detection 
by a spectrometer like SPI even more difficult. For 
that reason, detection of the ^^Ti lines from CasA 
and measurement of their width is one of the most 
critical issues in 7-ray line astronomy today, and one 
of the prime objectives of 5F//INTEGRAL. 

Finally, continuing search for ^''Ti lines in the Milky 
Way could reveal young supernova remnants and put 
interesting statistical constraints on ^"'Ti yields and 
supernova frequencies. Such a search has been al- 
ready performed in the past with IIEAO-3, SMM and 
COMPTEL data; a Monte Carlo analysis by The et 
al. (2000) showed that a rather rare supernova type 
with high '*'*Ti yield is favoured by the absence of 
a "^^Ti signal from the inner Galaxy. Similar conclu- 
sions are reached with preliminary INTEGRAL data 
(see contribution by Renaud et al. this volume). 



3. 26 AL AND 6"FE FROM MASSIVE STARS 

2^A1 is the first radioactive nucleus ever detected in 
the Galaxy through its characteristic gamma-ray line 
signature, at 1.8 MeV (Mahoney et al. 1984). Taking 




Figure 7. Map of Galactic ^^Al, after 9 years of 
COMPTEL operation (Pliischke et al. 2001a). 



into account that its lifetime of ^1 Myr is short w.r.t. 
galactic evolution timescales, its detection convinc- 
ingly demonstrates that nucleosynthesis is still active 
in the Milky Way (Clayton 1984). The detected flux 
is ~4 10"'' cm~2 s""' and corresponds to ~2 Mq of 
2^A1 currently present in the ISM (and produced per 
Myr, assuming a steady state situation). 

The COMPTEL instrument aboard CGRO mapped 
the 1.8 MeV emission in the Milky Way and found it 
to be irregular, with prominent "hot-spots" probably 
associated with the spiral arms (Diehl et al. 1995). 
The spatial distribution of ^^Al suggests that mas- 
sive stars are at its origin (Prantzos 1991, Prantzos 
and Diehl 1996), since any population of long-lived 
objects (AGB stars or novae) is expected to be more 
smoothly distributed. 

The detailed mapping of the Galactic distribution 
of 2^A1, obtainable through a determination of the 
distances to the "hot-spots" (see Kretsclimer et al. 
2003), is one of the main long-term objectives of IN- 
TEGRAL, since it will provide the most accurate 
picture of recent star formation in the Milky Way 
(see contributions by Diehl et al. and Hartmann et 
al., this volume). On the other hand, the study of in- 
dividual "hot-spots" revealed by COMPTEL bears 
on our understanding of the evolution of young stel- 
lar associations (in the cases of Cygnus, Carina and 
Centaurus-Circinus) and even individual stars (in the 
case of Vela). 

The Cygnus region was studied with population syn- 
thesis models by two groups (Cervinho et al. 2001, 
Pliischke et al. 2001b). The resulting morphology 
of the 1.8 MeV emission compares well with the 
COMPTEL data. However, in the case of Carina, 
the predicted absolute flux is smaller (by a factor 
of 5-20) than detected by COMPTEL (Knodlseder 
et al. 2001). That discrepancy may imply some- 
thing interesting, either for the (in)completeness of 
the stellar census of that association or for the ^^Al 
yields. INTEGRAL will establish more accurately 
the morphology of those "hot-spots" and further test 
the "massive star group" origin of 2^Al(see contribu- 
tion by Knodlseder et al. in this volume). 

The Vela region offers the opportunity to measure 
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Figure 8. Radial profiles of Al (solid curve) and 
Fe (dashed curve) in a 25 Mq star (from Rausher 
et al. 2002). The dotted curve is also for ^°Fe and 
it is obtained by reducing by a factor of 2 the cross 
section of Fe(n,j )^^ Fe (Al. Heger, private commu- 
nication). 

(or put upper limits on) ^^Al yields from individual 
sources. The morphology of the rather extended 1.8 
MeV emission detected by COMFTEL (Diehl 2002) 
does not allow identification with any of the three 
known objects in the field (the Vela SNR, the closest 
WR star 7^ Vel and SNR RX-J0852-4622); all three 
objects lie closer than 260 pc, according to recent es- 
timates. COMFTEL measurements are compatible 
with current yields of SNII (in the case of Vela SNR) 
and marginally compatible with current yields of 7 
Vel (Oberlack et al. 2000). INTEGRAL measure- 
ments in the Vela region are then expected to place 
more stringent constraints on stellar nucleosynthesis 
models (see contribution by Mowlavi et al. in this 
volume). 

It is not yet clear whether the majority of observed 
^^Al originates from the winds of the most massive 
stars (i.e. above 30 Mq, evolving as Wolf-Rayet 
stars) or from the explosions of less massive stars 
(i.e. in the 10-30 M© range, exploding as SNII); the 
uncertainties in the corresponding stellar yields are 
still quite large and do not allow to conclude yet. 

Clayton (1982) pointed out that SNII explosions 
produce another relatively short lived radioactivity, 
6°Fe (lifetime --2 Myr). Since WR winds do not 
eject that isotope, the detection of its characteristic 
gamma-ray lines in the Milky Way would constitute 
a strong argument for SNII being at the origin of 
^^Al. Based on detailed nucleosynthesis calculations 
of SNII (from Woosley and Weaver 1995, hereafter 
WW95) Timmes et al. (1995) found that the ex- 
pected gamma-ray line flux ratio of ^'^Fe/^^Al (for 
each of the two lines of ^°Fe) is 0.16, if SNII are the 
only sources of ^^Al in the Milky Way. 

The Reuven Ramaty High Energy Solar Spectro- 
scopic Imager (RHESSI) detected the galactic ^^Al 
emission at a flux level compatible with previous ob- 



Figure 9. The expected ratio of^^Fe/^^Al decays (for 
each of the two ^^Fe lines), convolved with a Salpeter 
stellar Initial Mass Function, is shown as a function 
of the upper stellar mass limit of the convolution in- 
tegral. The four curves correspond to the four differ- 
ent sets of stellar yields, with their thick portions cor- 
responding to the mass range covered in those works 
(see text). The dotted horizontal line at 0.10 is the 
^°Fe/^^Al line flux ratio reported by RHESSI (Smith, 
this volume) or, perhaps, an upper limit on that ra- 
tio. Filled pentagons mark the upper mass in each of 
the four studies; all recent calculations predict much 
higher values than the older calculations of WW95 
or the value of RHESSI. Only by taking into account 
the '^^ Al yields of massive Wolf-Rayet stars (thin por- 
tion of the curves beyond the masses indicated by the 
filled pentagons, obtained by adding data on Alfrom 
WR stars from Meynet et al. 1997) one may obtain 
^^Fe/^^ Al ratios compatible with the RHESSI results. 



servations (Smith 2003). Recently, Smith (2004) re- 
ported the first ever detection of the Galactic ^'^Fe 
gamma-ray lines with RHESSI. A reanalysis of the 
data (see contribution by Smith, this volume) leads 
to a marginal significance of the detection (2.6 cr) and 
to a hue flux ratio ^^Fef^^Al of 0.10 (for each ^'^Fe 
line), i.e. slightly lower than predicted by Timmes 
et al (1995) on the basis of WW95 nucleosynthesis 
calculations. 

However, more recent studies of SNII nucleosyn- 
thesis (Rauscher et al. 2002, hereafter RHHW02; 
and Limongi and Chieffi 2003, hereafter LC03) pro- 
duce different yields than WW95 (more ^°Fe and less 
^^Al); in particular, a large amount of ^"Fe is pro- 
duced in the He-layer (Fig. 8), probably due to the 
different nuclear reaction rates used in the new calcu- 
lations. As a result, the ^"Fe/^^Al ratio turns out to 
be considerably higher than the one found in Timmes 
et al. (1995). 

Prantzos (2004) compared the new yields with obser- 
vations, after convolving with a stellar Initial Mass 
Function (IMF). The results are plotted in Fig. 9. 
In the case of LC03 two curves are shown: LC03H 
corresponds to the high ^°Fe yields (low explosion 
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Figure 10. '^^ Al yields of massive stars, either ejected 
by SNII explosions (from Rauscher et al. 2002 and 
Limongi and Chieffi, 2003) or expelled from WR 
winds (from Meynet et al. 1997). SNII yields are 
for stars of metallicity Zq, while those of WR for 
metallicities Z^and 2 Zq; the latter is more appro- 
priate for the average metallicity of the Milky Way. 



energies) and LC03L to the low ^°Fe yields (high ex- 
plosion energies). In all cases the thick portions of 
the curves correspond to the stellar mass range cov- 
ered by each study, while the filled hexagons mark 
the highest mass of each calculation and thus pro- 
vide the IMF weighted value over the whole mass 
range covered by each study. It is clear from Fig. 9 
that recent calculations are in disagreement with the 
RHESSI result. 

Assuming that the recent theoretical results are not 
to be substantially revised in the future and that 
the RHESSI result is confirmed, what are the implica- 
tions for our understanding of the origin of ^^Al? The 
obvious conclusion is that the bulk of galactic ^^Al, 
detected by various instruments including RHESSI, 
is not produced by the source of ^°Fe: if this were 
the case, then ^°Fe would be detected with a line 
flux similar to the one of ^^Al. Obviously, another 
source of ^^Al is required, producing much smaller 
60pg/26Al ratios than the SNII. 

The obvious candidate source is Wolf-Rayet stars, as 
has been argued in many places over the years (e.g. 
Dearborn and Blake 1985, Prantzos and Casse 1986, 
Prantzos 1993, Prantzos and Dichl 1996, Meynet et 
al. 1997, Knodlseder 1999). The winds of those mas- 
sive, mass losing stars, eject large amounts of ^^Al 
produced through H-burning in the former convec- 
tive core, before its radioactive decay. In stars with 
no mass loss, those quantities of Al decay inside the 
stellar core before the final explosion and never get 
out of the star. Unfortunately, no consistent model 
for the complete evolution of a massive, mass losing, 
star (i.e. up to the final explosion) exists up to now. 

Adopting the ^^Al yields of non-rotating WR stars 
of solar initial metallicity by Meynet et al (1997), 



Figure 11. Radial distributions of A126, star for- 
mation rate (SFR) and metallicity (Z) in the Milky 
Way disk. Upper panel: Data points with vertical 
error bars correspond to various tracers of the SFR, 
while the galactic metallicity profile of oxygen ( with 
a gradient of dlog(O/II)=-0.07 dex/kpc) is shown by 
a solid line; the A126 profile, after an analysis of 
COMFTEL data by Knodlseder (1997), is shown 
(in relative units) by data points with vertical and 
horizontal error bars (the horizontal ones correspond 
to the adopted radial binning). Lower panel: If galac- 
tic A126 originates mostly from WR stars, its ra- 
dial distribution should scale with SFR * Z^'^ (points 
with vertical error bars, scaled from the upper panel), 
since the A126 yields of WR stars scale with Z^'^ at 
least (Vuissoz et al. 2003, calculations for rotating 
stars); however, the observed A126 distribution (same 
points as in upper panel) is flatter than the expected 
one in that case. 

which concern stars in the 25-120 M0 mass range, 
and combining them with the aforementioned SNII 
yields (see Fig. 10), one obtains the ™Fe/^^Al hne 
flux ratio produced by the total mass range of mas- 
sive stars, during all the stages of their evolution; 
this is expressed in Fig. 9 by the continuation of the 
four theoretical curves above the masses indicated by 
the filled pentagons. It can be seen that the RHESSI 
result is recovered in that case, provided that at least 
half of ^^Al originates from WR stars (in the case of 
LC03L), or even that 80% of ^'^Al originates from 
WR stars (in the case of RHHW02 or LC03H). 

Those arguments point towards WR stars as major 
sources of ■^^Al in the Milky Way. However, the sit- 
uation is far from clear yet, because the WR stellar 
yields of ^^Al depend strongly on metallicity. In the 
case of non rotating stellar models that dependence 
is cx Z^, according to Meynet et al. (1997). The 
rotating models of WR stars, currently calculated 
by the Geneva group (Meynet and Maeder 2003) 
show that rotation considerably alleviates the need 
for high mass loss rates, while at the same time lead- 
ing to the production of even larger ^^Al yields than 
the non- rotating models (Vuissoz et al. 2003, Pala- 
cios et al. this volume); in that case, it is found that 
the ^^Al yields of WR have a milder dependence on 
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metallicity (oc Z^-^) than the non rotating ones. 

In both cases, that mctaUicity dependence of the ^^Al 
yields of WR stars, combined with the radial profiles 
of star formation rate (SFR) and of metallicity in the 
Milky Way (see Fig. 11, upper panel) suggest that 
the resulting radial profile of ^^Al should be much 
steeper than the one actually observed. The latter, 
derived from COMPTEL observations (Knodlseder 
1997) appears in Fig. 11 (lower panel) and is clearly 
flatter than the product SFR*Z'^ ''' (as already no- 
ticed in Prantzos 2002). Similar conclusions are 
reached if the longitude, rather than radial, profiles 
of ^^Al, metallicity and SFR are considered. 

Thus, almost twenty years after its discovery, the 
^^Al emission of the Milky Way has not yet found a 
completely satisfactory explanation. Indeed, the re- 
cent observational {COMPTEL, RHESSI) and theo- 
retical (RHHW02, LC03, Vuissoz et al. 2003) results 
have made the puzzle even more complex than be- 
fore. The solution will obviously require progress in 
both theory and observations. From the theory point 
of view, detailed nucleosynthesis calculations of mass 
losing and rotating stars up to the final explosion in 
the mass range 12-100 M0 and for metallicitics up 
to 3 Zq will be required ; furthermore, the uncer- 
tainties still affecting the reaction rates of ^^Ne(a,n) 
(major neutron producer during He burning in mas- 
sive stars) and ^^Fe(n,7)^°Fe will have to be substan- 
tially reduced. From the observational point of view, 
the radial distributions of both ^^Al and ^°Fe will be 
needed; such distributions will probably be available 
if the operation of INTEGRAL is prolonged for a few 
years beyond its nominal 2-year mission (see contri- 
bution by C. Winkler, this volume). 



4. POSITRON ANNIHILATION IN THE 
GALAXY 

The first 7-ray line ever detected outside the solar 
system was the 511 keV line of electron-positron an- 
nihilation (Johnson et al. 1972, Leventhal et al. 
1978). Observations by various instruments in the 
90 's established that the line is not variable (at least 
in a ~10 year period), that its spatial distribution is 
apparently dominated by a bulge-like component and 
that the overall spectrum suggests a large positro- 
nium fraction of 0.93 (see Kinzer et al. 2001 and 
references therein). The 511 keV flux detected by 7- 
ray spectrometers (HEAO-3, ORIS, HEXAGONE, 
TGRS) in the central Galactic sterad was found to 
be ~10^'^ ph/cm-^/s, corresponding to a steady state 
production rate of 10*^ e+ s~^. 

Observations during the early mission phase of IN- 
TEGRAL broadly confirm that picture. Analysis of 
SPI data (Knodlseder et al. 2003, Jean et al.' 2003, 
contribution by Jean et al. this volume) finds a 
rather narrow line (FWHM~3 keV), excludes a point 
source in the Galactic Center (GC), and finds that 
the emission is bulge-like (best fit with a Gaussian of 
FWHM=9°). The non-detection of a disk emission 



up to now imposes a model-dependent lower limit on 
the bulge/disk ratio of 0.4-0.8 (see also contribution 
by Weidenpointer et al. this volume). The detected 
flux of 9.9l2!i 10~^ ph/cm^/s is compatible with pre- 
vious measurements. 

Possible sources of Galactic positrons have been 
studied in several works up to now (e.g. Dermer and 
Murphy 2001 and references therein). Most promi- 
nent among them seems to be the -radioactivity of 
supernovae (see Table 1) and in particular of SNIa. 
On the other hand, the already detected emission of 
^^Al in the Galaxy (Sec. 3) implies that this radioac- 
tive nucleus contributes a non-negligible fraction of 
the Galactic positron production rate. A thorough 
analysis of the supernova contribution to positron 
production in the Galaxy has been done by Milne 
et al. (1999a). Their conclusion was that the decay 
of ■''''Co from the total Galactic population of SNIa 
may produce ^ 10*"^ e+ s^^, and that this rate could 
explain about half of the the total Galactic 511 keV 
emission as measured by OSSE, SMM and TGRS. 

The improved spatial (and spectral) resolution of SPI 
gives a dramatic new twist to the issue of Galactic 
positrons, because it is clear now that ^10^'^ e+ s^-"^ 
are required to be produced in the bulge alone, not in 
the whole Milky Way (unless a mechanism is found 
to channel all positrons from the disk to the bulge). 

The positron production rate from SNIa in the bulge 
is 

R = M F N 

where: M is the mass of the bulge, F is the SNIa fre- 
quency per unit mass in a bulge-like system, N is the 
number of positrons ejected by a typical SNIa. Typi- 
cal values for those parameters are: M=1.5to;5 ^0^° 
Mq (e.g. Launhardt et al. 2002), J^=0.22t^;g^ SNIa 
per 10^" M0 per millenium (Cappellaro et al. 2003, 
assuming that the bulge is of intermediate type, be- 
tween EO and Sa) and N=?,t\ 10^^ positrons (Milne 
et al. 1999b, i.e. the ^^Co positron escape fraction 
is ~3%). With those numbers one gets R ^ 0.8 10^^ 
e"*" s""'^, i.e. a production rate about 12 times lower 
than suggested by the SPI data analysis. The ex- 
pected contribution from other radioactivities (^''Ti 
from SN, ^^Al from massive stars, ^^Na from no- 
vae) is too low to account for the missing amount of 
positrons. 

It should be stressed that the uncertainties on F and 
N could be considerably larger than quoted in the 

previous paragraph. Indeed, the fact that SNIa have 
not been detected up to now in the bulges of spirals 
(Wang et al. 1997) suggests that important selec- 
tion biases, e.g. extinction, may affect the results of 
current surveys. 

On the other hand, the value of N has been em- 
pirically derived from the late lightcurves of SNIa 
(Milne et al. 1999b). This estimate is on much 
more secure grounds than the theoretical estimates of 

Chan and Lingelfclter (1993). However, it may also 
be affected by uncertainties on any early escape of 
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Table 2. Candidate sources for positrons in Galactic bulge 



SOURCE 


MORPHOLOGY 


INTENSITY 


COMMENTS 


SNIa 
radioactivity 


(,n'nQiQ"t"A'n1" ^s^il"}! 

v^UiiOio LV^iiU WlLrii 

observations 


( "standard" parameter values) 


1 1 n PPT"t" fll Tli"! PQ 

(perhaps) underestimated 


LMXB 
Outflows / ip1"s 


Seems OK, BUT 

nrnninjitprl nv tpw ^trnnp' 

sources outside the bulge 


Unknown, but 

nop^ Tint ^ppTTi nnTP^mnnJi nip 

because — > 


Less than 1% of available 
to form positrons 


Hypernova(e) 


Improbable 
(how to fill the bulge 
with positrons?) 


Not unreasonable, 
but very poorly 
constrained 


More hypernovae 
expected in disk 
(molecular ring) 


Dark Matter 
annihilation 


Centered on GC, 
otherwise unknown 


Unknown/ 
Unconstrained 


Postulated light particle 
should be seen in LEP 
(if in Standard Model) 



pared to ~10'^^ erg/s required to produce ^10^'^ e"*" 
s^^ with energies of ^^1 MeV per positron, as those 
produced from radioactivity. It is sufficient to con- 
vert only ~1 % of the available energy into positrons 
to explain the bulge emission. However, the mech- 
anism of that conversion is not known yet. More- 
over, the observed X-ray emissivity is dominated by 
a dozen strong sources, which are mostly distributed 
in the Galactic plane and not in the bulge (see Table 
5 in Grimm et al. 2002). 

Even more speculative is the suggestion (Nomoto et 
al. 2001) that very energetic (and presumably as- 
pherical) explosions of massive stars, known as "hy- 
pernovae" are at the origin of the 511 keV emision. 
This idea has been further elaborated in Casse et 
al. (2004, see also contribution by Schanne et al. 
this volume). The observed constraints on 511 keV 
intensity could be satisfied in that case, since such 
explosions are expected to produce large amounts 
of ^^Ni (see Fig. 3, and SN1998bw in particular) 
along the rotation axis, thus making easier the es- 
cape of positrons. However, there is no quantita- 
tive evaluation yet of the positron escape fraction 
in such sources, and any theoretical one will suf- 
fer from the same large uncertainties as the analy- 
sis of Chan and Lingenfelter (1993) on SNIa; em- 
pirical evaluation will require much larger statistical 
samples and longer observations of hypernovae than 
presently available. Independently of that, the idea 
has several shortcomings: it is difficult to imagine 
how the positrons of a single (or a few) explosion(s) 
could fill the bulge; and it is statistically improba- 
ble that such events are not detected in the Milky 
Way disk (for instance in the molecular ring), where 
the star formation rate is much higher than in the 
galactic center or bulge and where all conditions for 
positron slow down and annihilation are fulfilled. 

Finally, annihilation of a rather special kind of light 
dark matter particles has been recently proposed as 
the source of galactic positrons (Boehm et al. 2003, 
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Figure 12. Longitude distribution of galactic LMXBs 
(thick histogram) and HMXBs (thin histrogram) , 
from Grimm et al. (2002). The collective emissiv- 
ity of LMXBs is ^10 times higher than the one of 
HMXBs and 100 times higher than required to pro- 
duce 10^'^ e"*" s~^. 



positrons from ^^Ni mixed in the outer layers. Thus, 
SNIa should not yet be definitively excluded as major 
sources of bulge positrons. However, other positron 
sources should be sought. 

Among the various candidate sources. Low Mass 
X-Ray Binaries (LMXBs) appear quite promising, 
at least as far as the morphology and the inten- 
sity constraints of the observed 511 keV emission 
are concerned. Indeed, the observed distribution of 
LMXBs in the Milky Way is strongly concentrated 
towards the Galactic bulge (Fig. 12), and satis- 
fies the bulge/disk ratio obtained by the prelimi- 
nary SPL analysis (P. Jean, private communication). 
Moreover, the collective X-ray emissivity of Galactic 
LMXBs is -10^^ erg/s (Grimm et al. 2002), com- 
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contribution by Casse et al. this volume). The pro- 
posed particles are quite light (in the 1-100 MeV 
range) so that their annihilation does not produce 
undesirable high energy gamma-rays, and in that 
respect they do not correspond to the most com- 
monly discussed dark matter candidate, which has 
mass in the GeV to TeV range. Moreover, rather 
special properties are required for such light parti- 
cles to justify why they have escaped detection up 
to now in accelerators such as the LEP. It is hard 
to evaluate the plausibility of that hypothesis, since 
the required properties of the source (i.e. density 
profile, annihilation cross-section) are completely un- 
known/unconstrained; in fact, the observed proper- 
ties of the 511 keV emission (intensity and density 
profile) are used in Boehm et al. (2003) in order to 
derive the properties of the dark matter source of 
positrons. 

Further observations by INTEGRAL will help, since 
it is expected that the Galactic disk will sooner or 
later manifest itself in the light of 511 keV photons. 
Premiminary hints for disk 511 keV emission were 
found in OSSE's data (e.g. Kinzer et al. 2001, 
Milne et al. 2002). The observed disk emission of 
^^Al at 1.8 MeV puts a lower limit on the positron 
production rate of the disk; observations will estab- 
lish whether those positrons annihilate locally (i.e. 
whether the steady-state assumption is locally valid) 
or whether they escape in the halo to annihilate at 
higher latitudes. In that respect, the Cygnus region, 
one of the most prominent "hot-spots" in COMP- 
TEL's 1.8 MeV map (see Fig. 7) is a most promising 
target. 
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